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ABSTRACT 

Eclipse lightcurves of the dwarf nova IP Peg during the November 1996 outburst 
are analysed with eclipse mapping techniques to constrain the location and investigate 
the spatial structure of the spiral shocks observed in the Doppler tomograms (Har- 
laftis et al. 1999). Eclipse maps in the blue continuum and in the CIII+NIII A4650 
emission line show two asymmetric arcs of ~ 90 degrees in azimuth and extending 
from intermediate to the outer disc regions (i? ~ 0.2 — 0.6 Rli, where R^i is the dis- 
tance from disc centre to the inner Lagrangian point) which are interpreted as being 
the spiral shocks seen in the Doppler tomograms. The Hell A4686 eclipse map also 
shows two asymmetric arcs diluted by a central brightness source. The central source 
probably corresponds to the low-velocity component seen in the Doppler tomogram 
and is understood in terms of gas outflow in a wind emanating from the inner parts 
of the disc. We estimate that the spirals contribute about 16 and 30 per cent of the 
total line flux, respectively, for the Hell and CIII-t-NHI lines. Comparison between 
the Doppler and eclipse maps reveal that the Keplerian velocities derived from the 
radial position of the shocks are systematically larger than those inferred from the 
Doppler tomography indicating that the gas in the spiral shocks has sub-Keplerian 
velocities. We undertake simulations with the aim to investigate the effect of artifacts 
on the image reconstruction of the spiral structures. 
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1 INTRODUCTION 

Accretion discs are widespread in astrophysical environ- 
ments, from sheltering the birth of stars to providing the 
energetics for the most violent phenomena such as relativis- 
tic jets. Despite its general importance and although con- 
siderable effort both in observation and theory has been in- 
vested over the past decade, the structure and underlying 
physics of accretion discs remains poorly understood. One 
of the major unsolved problems concerns the nature of the 
anomalous viscosity mechanism responsible for the inward 
spiraling of the disc material (Frank, King & Raine 1992). 

Best prospects for progress in understanding accretion 
discs physics are possibly found in mass-exchanging binaries 
such as Cataclysmic Variables (CVs). In these close binaries 
mass is fed to a non-magnetic {B ^ 10"" G) white dwarf via 
an accretion disc by a Roche lobe filling companion star 



(the secondary). The sub-class of dwarf novae comprises 
low-mass transfer CVs which show recurrent outbursts of 
2-5 magnitudes on timescales of months either due to an 
instability in the mass transfer from the secondary or due 
to a thermal instability in the accretion disc which switches 
the disc from a low to a high-viscosity regime (Warner 1995 
and references therein). 

Spiral shocks have been advocated by various re- 
searchers as a possible mechanism for transport of angular 
momentum in accretion discs (Savonije, Papaloizou & Lin 
1994) and may be the key, together with magnetic viscos- 
ity (Hawley, Balbus & Winters, 1999) , in understanding the 
viscosity mechanism. The recent discovery of spiral shocks 
in the accretion disc of the dwarf novae IP Pegasi in out- 
burst - from Doppler tomography of emission lines (Steeghs, 
Harlaftis & Home 1997, 1998; Harlaftis et al. 1999) - con- 
firmed the results of hydrodynamical simulations (Armitage 
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& Murray 1998, Stehle 1999). The spiral shocks are pro- 
duced in the outer regions of the disc by the tides raised by 
the secondary star. During the outburst the disc expands 
and its outer parts feel more effectively the gravitational at- 
traction of the secondary star leading to the formation of 
spiral arms. 

Here we report on the eclipse mapping analysis of the 
data obtained by Harlaftis et al. (1999; see there for ob- 
servations and data reduction). Our goal is to confirm the 
existence, constrain the location and to investigate the spa- 
tial structure of the spiral shocks observed in the Doppler 
tomograms. Section^ presents the data and gives details of 
the analysis procedures. In section^ we present a set of sim- 
ulations with the eclipse mapping method aimed to clarify 
the interpretation of the results of section^ in terms of real 
spiral shocks. A summary of our findings is given in section^ 



2 DATA ANALYSIS 
2.1 Lightcurves 

A time-series of high-resolution, optical spectrophotometry 
(AA = 4354 - 4747 A, velocity dispersion of 27 fcms"^ per 
pixel) covering one eclipse of IP Peg was obtained during 
the third day of the November 1996 outburst. The reader is 
referred to Harlaftis et al. (1999) for a detailed description 
of the dataset and of the reduction procedures. Lightcurves 
were extracted for the blue continuum (4365 — 4440 A) and 
for the CIII-HNIII A4650 (Bowen blend) and Hell A4686 
lines and phase folded according to the sinusoidal ephemeris 
of Wolf et al. (1993), 

rmid(HJD) = 2 445 615.4156 + 0.158 206 16 £ + (0-C) (1) 



where 



(O - C) = 1.0903 X 10"^ sin 



2tt 



E - 10258 



10850.88 



The line lightcurves were continuum subtracted and, there- 
fore, correspond to net line emission. The three lightcurves 
are shown in Fig. ^ as gray open squares. 

The C III-I-N HI lightcurve shows a peculiar double- 
stepped eclipse shape revealing the presence of two asym- 
metric brightness sources displaced from disc centre. Al- 
though less pronounced, the same morphology can also be 
seen in the continuum lightcurve. The shape of the Hell 
eclipse is more symmetrical than that of the other passbands 
but mid-eclipse occurs earlier with respect to the continuum 
eclipse, indicating that the line surface distribution is also 
asymmetric. 

The continuum and Hell lightcurves show a conspicu- 
ous orbital modulation with maximum at phase (j> ~ —0.15 
cycle and minimum at ~ -1-0.15 cycle. This is not seen 
in the C III-I-N HI lightcurve although an increase in flux is 
clearly visible after phase (j) = -1-0.22 cycle. For any reason- 
able mass ratio, g < 1 [(?=0.5 for IP Peg, Wood & Craw- 
ford (1986)], the shadow of the secondary star covers re- 
gions outside the primary lobe for orbital phases |0| > 0.2 
cycle. Therefore, it is hard to explain the observed modu- 
lation in terms of occultation by the secondary star unless 
the eclipsed source lies outside the primary lobe. Thereafter, 
we assign the orbital modulation to gas obscuration by the 



spiral arm seen at maximum between phases 0.0-0.25 cycle 
(see section^). 

Out-of-eclipse brightness changes are not accounted for 
by the eclipse mapping method, which assumes that all vari- 
ations in the eclipse lightcurve are due to the changing occul- 
tation of the emitting region by the secondary star (but see 
Bobinger et al. 1996 for an example of how to include orbital 
modulations in the eclipse mapping scheme). Orbital varia- 
tions were therefore removed from the lightcurves by fitting 
a spline function to the phases outside eclipse, dividing the 
lightcurve by the fitted spline, and scaling the result to the 
spline function value at phase zero. This procedure removes 
orbital modulations with only minor effects on the eclipse 
shape itself. The corrected lightcurves are shown in Fig. nl as 
filled circles. For the purpose of eclipse mapping analysis the 
lightcurves were limited to the phase range (—0.18, -1-0.28) 
since data outside of eclipse is basically used to set the out- 
of-eclipse level and fiickering in this phase range serves to 
add unnecessary noise to the fit. 

2.2 Eclipse maps 

The eclipse mapping method is an inversion technique that 
uses the information contained in the shape of the eclipse to 
recover the surface brightness distribution of the eclipsed ac- 
cretion disc. The reader is referred to Horne (1985), Rutten 
et al. (1992) and Baptista & Steiner (1993) for the details 
of the method. 

As our eclipse map we adopted a grid of 51 x 51 pix- 
els centered on the primary star with side 2 Rli, where 
Rli is the distance from the disk center to the inner La- 
grangian point. The eclipse geometry is controled by a pair 
of q and i values. The mass ratio q defines the shape and 
the relative size of the Roche lobes. The inclination i deter- 
mines the shape and extent of the shadow of the secondary 
star as projected onto the orbital plane. We obtained recon- 
structions for two sets of parameter, (g = 0.5 , i — 81°) 
(Wood & Crawford 1986) and (g = 0.58 , i = 79.5°) (Marsh 
1988), which correspond to an eclipse width of the disc 
centre of A0 = 0.086 (Wood & Crawford 1986; Marsh & 
Horne 1990). These combination of parameters ensure that 
the white dwarf is at the center of the map. There is no 
perceptible difference in eclipse maps obtained with either 
geometry. Hence, for the remainder of the paper we will refer 
to and show the results for (g — 0.5 , i — 81°). 

The lightcurves were analyzed with eclipse mapping 
techniques to solve for a map of the disc brightness distribu- 
tion and for the ffux of an additional uneclipsed component 
in each passband. The uneclipsed component accounts for all 
light that is not contained in the eclipse map in the orbital 
plane (i.e., light from the secondary star and/or a vertically 
extended disc wind). The reader is referred to Rutten et al. 
(1992) and Baptista, Steiner & Horne (1996) for a detailed 
description of and tests with the uneclipsed component. For 
the reconstructions we adopted the default of limited az- 
imuthal smearing of Rutten et al. (1992), which is better 
suited for recovering asymmetric structures than the orig- 
inal default of full azimuthal smearing (see Baptista et al. 
1996). 

Lightcurves, fitted models and grayscale plots of the re- 
sulting eclipse maps are shown in Fig. Hand will be discussed 
in detail in sectionM. 
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3 ECLIPSE MAPPING SIMULATIONS 

We performed various simulations with asymmetric sources 
in order (i) to investigate how the presence of spiral struc- 
tures in the accretion disc affects the shape of the eclipse 
lightcurve, and (ii) to evaluate the ability of the eclipse map- 
ping method to reconstruct these structures in the eclipse 
maps. 

For the simulations we adopted the geometry of IP Peg 
(q = 0.5 and i = 81°) and constructed lightcurves with the 
same signal to noise ratio and orbital phases of the real data 
of section p. Figure^ shows the results of the simulations. 

Asymmetric compact sources (model #1) result in 
eclipse lightcurves with rapid brightness changes at 
ingress/egress phases. The azimuthal smearing effect charac- 
teristic of the eclipse mapping method is responsible for the 
distortion which makes the compact sources appear 'blurred' 
in azimuth. Nevertheless, their radial and azimuthal loca- 
tions are satisfactory recovered. 

Brightness distributions with spiral structures (mod- 
els #2 to #4) result in eclipse shapes with characteristic 
'bulges', whose extension and location in phase reflect the 
orientation and radial extent of the spiral arms. Due to the 
azimuthal smearing effect these structures are reproduced 
in the form of asymmetric arcs, whose maximum brightness 
and radial position yield information about the orientation, 
position and radial extent of the original spiral arms. The 
adition of a symmetric brightness source (i.e., centred in the 
eclipse map, model #5) can dilute the presence of spiral 
arms. In this case the eclipse shape is smoother and more 
symmetric in comparison with those of models #2-4 and the 
asymmetric arcs are less clearly visible in the eclipse map, 
mixing with the brightness distribution of the central source. 

These simulations show that the eclipse mapping 
method is able to reproduce asymmetric light sources such 
as spiral arms (provided the asymmetric sources are prop- 
erly eclipsed) and that the asymmetric structures seen in 
the eclipse maps of Fig.|^ are not caused by artifacts of the 
method. Models #3 to #5 are the relevant ones for the pur- 
pose of comparing the results of the simulations with those 
from the IP Peg data (Fig.|^. The morphology of the con- 
tinuum and C III-I-N III lightcurves is similar to that of mod- 
els #3 and #4, while the Hell lightcurve resembles that of 
model #5. 



4 RESULTS 

Data and model lightcurves are shown in the left panels of 
Fig.|2|. Horizontal dashed lines indicate the uneclipsed com- 
ponent in each case. The uneclipsed component corresponds 
to about 12, 1 and 1 per cent of the total flux, respectively 
for the continuum, C III-I-N HI and Hell curves. While a 
non-negligible fraction of the light in the continuum prob- 
ably arises from an emitting region outside of the orbital 
plane (possibly a disk wind), the net Hell and C III-I-N HI 
emission mostly arises from (or close to) the orbital plane. 

The middle panels of Fig. ^ show eclipse maps in a 
logarithmic grayscale. The maps in the right panels show 
the asymmetric part of the maps in the middle panels and 
are obtained by calculating and subtracting azimuthally- 
averaged intensities at each radius. 



The continuum and C III-I-N HI lightcurves display 
bulges similar to that of the lightcurves of models #3 and 
#4 (Section 3) and result in eclipse maps with two clearly 
visible asymmetric arcs, which are interpreted as being the 
spiral shocks seen in the Doppler tomograms of Harlaftis 
et al. (1999). In comparison with the models of Fig. ^, the 
orientation of the arcs suggests that the spirals are aligned 
in a direction perpendicular to the major axis of the binary 
(models ^^3 and #4). The arcs show an azimuthal extent of 
~ 90° and extend from the intermediate to the outer disc 
regions {R ~ 0.2 — 0.6 -Rli). Therefore, the outer radius of 
the spirals is of the same order of the disc outburst radius 
(Rd ~ 0.34 a ~ 0.6 Rli) inferred by Wood et al. (1989). The 
eclipse maps show no evidence of the bright spot at disc rim 
and no enhanced emission along the gas stream trajectory. 

The azimuthal location of the arcs is consistent with 
the results from hydrodynamical simulations and from the 
Doppler tomography. The arc in the upper left quadrant of 
the eclipse map (hereafter arc 1) corresponds to the spiral 
arm whose maximum occurs at phases 0.5-0.75 while the arc 
in the lower right quadrant (arc 2) corresponds to the spiral 
arm seen at maximum intensity at phases 0.0-0.25 (orbital 
phases increases clockwise in the eclipse maps of Fig. ^ and 
phase zero coincides with the inner lagrangian point LI). 
The arcs are not symmetrical with respect to the centre of 
the disc. In C III-I-N HI, arc 2 is further away from disc centre 
than arc 1 - in agreement with the Doppler tomography, 
which indicates smaller velocities for the spiral arm 2 (~ 
550 kms~^) than for the spiral arm 1 (~ 700 kms~^). 

The lightcurve of Hell is quite symmetrical with less 
pronounced bulges than in C III-I-N HI, resulting in an 
eclipse map consisting of a symmetrical, centred brightness 
distribution and asymmetric arcs at different distances from 
disc centre. The outermost arc (arc 2) is more easily seen in 
the eclipse map, while the emission from the innermost arc 
(arc 1) is blendend with that of the central source. Never- 
theless, arc 1 is clearly seen in the asymmetric part of the 
He II map show in the right panel of Fig.|2| The symmetrical 
emission component is probably related to the low-velocity 
component seen in Hell Doppler tomograms and is sug- 
gested to be due to gas outflow in a wind emanating from 
the inner parts of the disc (see also Marsh and Horne, 1990; 
for an alternative interpretation, slingshot prominence from 
the secondary star, see Steeghs et al. 1996). The Hell arcs 
contribute about 16 per cent of the total flux of the eclipse 
map - in good agreement with the results from the Doppler 
tomography, which indicate that the spirals contribute ~ 15 
per cent of the total Hell emission (Harlaftis et al. 1999). 
In comparison, the arcs in the CIII4-NIII and continuum 
maps contribute, respectively, about 30 and 13 per cent of 
the total flux. 

We quantify the properties of the asymmetric arcs by 
dividing the eclipse map in azimuthal slices (i.e., 'slices of 
pizza') and computing the distance at which the intensity 
is maximum for each azimuth. This exercice allows to trace 
the distribution in radius and azimuth of the spiral struc- 
tures. The results are plotted in Fig.^ as a function of or- 
bital phase. The diagrams for Hell were computed from its 
asymmetric map and are noisier than those for C III-I-N HI 
because most (~ 84 per cent) of the flux in the eclipse map 
is subtracted with the symmetric component. The two spi- 
ral shocks are clearly visible in the intensity diagrams, as 



© 2000 RAS, MNRAS 000, 1-7 



4 R. Baptista et al. 



well as their distinct locations with respect to disc centre. 
In He II the outer spiral (arc2) is brighter than the inner spi- 
ral (arcl), in line with the results of Harlaftis et al. (1999), 
while in CIII+NIII arc 1 is brighter than arc 2. The middle 
panels give the radial position of the maximum intensity as a 
function of binary phase. For CIII+NIII, the maximum in- 
tensity along arc 1 lies at a constant distance of ~ 0.28 Rli 
from disc centre while the maximum intensity of arc 2 occurs 
at ~ 0.55 Rl\ - The numbers are similar for Hell. 

We computed equivalent Keplerian velocities for each 
radius assuming Mi = 1.0 ± 0.1 Mq e i?Li = 0.81 Rq 
(Marsh & Horne 1990). The results are plotted in the upper 
panel of Fig. ^. Gray lines show the corresponding uncertain- 
ties at the 1-a limit. For comparison, the results from the 
CIII-I-NIII and Hell Doppler tomograms (Harlaftis et al. 
1999; see their fig. 4 for the Hell diagram) are shown as 
dashed lines. Since the CIII-I-NIII Doppler map is much 
noisier and blurred than the Hell Doppler map, the corre- 
sponding diagram is noisier and less reliable than the He II 
diagram on the right panel. We obtain velocities in the range 
850 — 1050 kms~^ for the spiral 1 compared to the observed 
400 -770 km and in the range 650 - 800 km for the 
spiral 2 compared to the observed 400 — 550 kms~^ (ob- 
served values from Harlaftis et al. 1999). The Keplerian ve- 
locities calculated from the radial position of the shocks are 
systematically larger than those inferred from the Doppler 
tomography, suggesting that the gas in the spiral shocks has 
sub-Keplerian velocities. This is in line with the results of 
the hydrodynamical simulations of Steeghs & Stehle (1999, 
see their fig. 5), which predicts velocities lower than Keple- 
rian (by as much as 15 per cent) in the outer disc near the 
spirals. 

We remark that with the white dwarf mass and Roche 
lobe radius of IP Peg the Keplerian velocity at the largest 
possible disc radius {R ~ 0.85 -Rli) is about 530 kms~^. 
Therefore, if the observed velocities of Harlaftis et al. (1999) 
do refiect Keplerian motions, then the emitting gas should 
be at the border and even outside the primary lobe. 

Occultation of light from the inner disc regions by the 
spirals might produce the out of eclipse variations seen in 
Fig.|l|. This is expected since the spiral waves are also ver- 
tically extended. From the azimuthal position of the spirals 
in the eclipse map, the maximum occultation (i.e., the min- 
imum of the orbital modulation) should occur when the spi- 
rals are seen face on, at orbital phases ~ —0.3 and ~ +0.15 
cycle, in agreement with the modulations seen in Fig.nl. 



5 CONCLUSIONS 

We analyzed eclipse lightcurves of the dwarf novae IP Peg 
during the November 1996 outburst in order to confirm the 
existence, constrain the location and investigate the spatial 
structure of the spiral shocks observed in the Doppler tomo- 
grams. Our mais results can be summarized as follows: 

• Eclipse maps in the blue continuum and in the 
C III+N HI emission line reveal two asymmetric arcs at dif- 
ferent azimuth and radius from disc centre which are con- 
sistent with the spiral shocks seen in the Doppler tomo- 
grams. The arcs show an azimuthal extent of ~ 90° and 
extend from the intermediate to the outer disc regions {R ~ 



0.2 — 0.6 Rli)- The outer radius of the spirals is of the same 
order of the disc outburst radius {Rd — 0.34 a ~ 0.6 Rli)- 

• The Hell eclipse map is composed of a central bright- 
ness source plus asymmetric arcs at different distances 
from disc centre. The symmetric component probably corre- 
sponds to the low- velocity component seen in He II Doppler 
tomograms and is understood in terms of gas outflow in a 
wind emanating from the inner parts of the disc. 

• The spirals contribute about 16 and 30 per cent of the 
total line flux, respectively, for the He II and C III+N III 
lines, and 13 per cent in the continuum. 

• The Keplerian velocities derived from the radial posi- 
tion of the shocks are systematically larger than those in- 
ferred from the Doppler tomography, indicating that the gas 
in the spiral shocks has sub-Keplerian velocities. 
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Figure 1. Original liglitcurves (gray open squares), fitted splines (gray dashed lines) and corrected lightcurves (black filled squares), 
for the continuum, Clll+Nin and He II data. Vertical dotted lines mark the ingress/egress phases of the white dwarf for an assumed 
eclipse width of = 0.086 cycle (Wood & Crawford 1986) and horizontal dotted lines mark the reference (mid-eclipse) flux level of the 
spline fit. 



Figure 2. Left: Data (dots with error bars) and model (solid lines) lightcurves of IP Peg at outburst maximum for the continuum, 
CIII+NIII A4650 and Hell A4686 lines. Horizontal dashed lines indicate the uneclipsed component in each case. Middle: eclipse maps 
in a logarithmic grayscale. Right: the eclipse maps of the middle panel after subtracting their symmetric part; these diagrams emphasize 
the asymmetric structures. Brighter regions are indicated in black; fainter regions in white. A cross mark the center of the disc; dotted 
lines show the Roche lobe and the gas stream trajectory; dotted circles mark disc radii oi R = 0.2 and 0.6 Rli] the secondary is to the 
right of each map and the stars rotate counterclockwise. 
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lightcurve model reconstruction 



#1 




-0.2-0.1 0.1 0.2 
orbital phase 

Figure 3. Simulations with asymmetric brightness distributions. Panels in the center show different synthetic images in logarithmic 
grayscale. Panels in the left show noise-added lightcurves derived from these brightness distributions (dots) and the fitted models (solid 
lines). Vertical dotted lines mark the ingress/egress phases of the center of the disc. The dashed curve in the lower panel illustrates the 
contribution from the spiral arms to the eclipse shape of model #5. Panels on the right show the corresponding eclipse maps in the same 
logarithmic grayscale as in the middle panel. The notation is the same as in Fig.bl 



Figure 4. The dependency with binary phase of the maximum intensity, radius and corresponding Keplerian velocity at maximum 
intensity, as derived from the CIII+NIII and Hell eclipse maps. The velocities were computed assuming Mi = 1.0 ± 0.1 Mq and 
R^yl = 0.81 i?0 (Marsh & Horne 1990) and the intensities are plotted in an arbitrary scale. Gray lines show the uncertainties in the 
velocity at the 1-cr limit. The dependency of velocity with orbital phase as derived from the CIII+NIII and Hell Doppler tomograms 
(Harlaftis et al. 1999; see their fig. 4) are shown as dashed lines for comparison. The location of the spiral arms are indicated by horizontal 
bars with labels 1 (inner spiral) and 2 (outer spiral). 
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